Dark Energy: Beyond General Relativity? 
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Abstract. The late-time accelerated expansion of the universe is a major challenge for cosmology. It may well be that 
a solution to this problem will require a theory of gravitation beyond General Relativity. It is emphasized that precision 
cosmology will strongly constrain the possibilities by using observational data probing the background as well as the 
inhomogeneities. 



INTRODUCTION 

There is growing observational evidence for a late time accelerated expansion of our universe d. Such an expansion 
constitutes a radical departure from concentional cosmology but it is reassuring that a consistent picture of our universe 
seems to emerge from all the data |2]. This accelerated expansion might be due to a new component with sufficiently 
negative pressure, coined Dark Energy. The simplest candidate of this kind is a positive cosmological constant A 
whose pressure satisfies p^ = — p*. As well-known, its interpretation as the vacuum energy is problemetic because of 
its exceeding smallness (see e.g. yl for recent comprehensive reviews). One can also introduce phenomenologically 
a smooth component with constant equation of state and sufficiently negative pressure typically satisfying woe = 
t^E- < —0.5 1 4]. Observations could force us to consider a smooth component with a varying equation of state. The 
most prominent candidate in this respect is a minimally coupled scalar field, Quintessence models [5]. A further 
constraint appears if observations force us to consider DE models of the phantom type, i.e. satisfying wde < — 1 l£j 
as observations seem to suggest on small redshifts z < 0.5 |7|]. The most striking consequence in that case is that 
quintessence models are ruled out. 

However, the late-time acceleration might as well be caused by a change in the theory of gravitation, no longer 
described by General Relativity |^|. Ineterest for such candidates is increased by their ability to produce DE of the 
phantom type. We will review scalar-tensor DE models 0, full as a promising candidate belonging to this class of 
models . It is emphasized that observations probing the background and the inhomogeneities will establish whether 
such attempts provide a consistent description of our universe. 



TWO BASIC DARK ENERGY MODELS INSIDE GENERAL RELATIVITY 
Constant equation of state DE models 

The simplest DE models are those in the framework of GR containing some unknown component with negative 
equation of state parameter wde = j^j ■ Our starting point are the well-known Friedmann equations for a homogeneous 
and isotropic universe, 
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In the equations above, the different components labelled i, are all isotropic perfect fluids while a dot stands for a 
derivative with respect to (cosmic) time t. We note from Q that a necessary, but not yet sufficient, condition that a 
component i induces an accelerated expansion is given by 



Pi + 3pi<0. (3) 

We concentrate on equations of state of the form p, = w, /?,-, with constant w, and we specialize to a universe containing 
dustlike matter and some DE component. From (|2}, it is not hard to deduce that our universe is presently accelerating 
provided 

w DE <-Ul + -^), (4) 
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and in particular for a flat universe 

WDE < -^D£,0 ■ ( 5 ) 

For instance, for i2,„,o = 0.3 , Q.de ,o = 0.7, wde < —0.47 is required. Therefore, present experimental evidence yielding 
^m,o ~ 0.3, combined with the location of the first acoustic peak in the CMB anisotropy first detected by balloon 
experiments and confirmed by the recent CMB data released by WMAP 1 12] suggesting a nearly flat universe, imply 
that our universe would be presently accelerating for a wide range of constant values, roughly woe < —0.5. Introducing 
the dimensionless quantity x = accelerated expansion starts at x a given by 

x a = (-1 +3|vv D£ |) — , (6) 
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which corresponds to redshifts 
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z a = (-l+3\w DE \)^DE\ -pi -1. (7) 
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For — 1 < wde — * —1/3, Z a is shifted towards smaller redshifts, for (Cl,„,o, Qde,o) — (0.3, 0.7), we have z a = 0.671 
for a constant A-term, and z a = 0.414 when wde = —0.6. The fact that z a is so close to zero, is the cosmic coincidence 
problem. Finally, these models can yield a significant increase, depending on WD£,Qm,o,QD£,o> of the age of the 
universe for given Hubble parameter Hq compared to an Einstein-de Sitter universe, i.e. a flat universe with £l m q = 1 



Quintessence 

The DE component could well be a time dependent minimally coupled scalar field 0(f) called Quintessence. This 
possibility is clearly inspired by the inflatonary paradigm in which a scalar field is so successful in implementing the 
inflationary stage. Such a scalar field can be considered as a perfect fluid with 

P^ = \f + V($) P^ = \^-V(^) (8) 

and therefore the equation of state parameter w a is given by 

w * = (9) 

For Ofj, > 0, the equation of state must satisfy 

vt> > -1 , (10) 

in other words (j) cannot be of the Phantom type. It is possible to have scaling solutions with x m , m = 

— 3(1 + w ( | ) )=constant. However this requires a very particular potential V (<p) for which, 

vW = lz^fl. (11) 

yv > l+w 2 V 



Hence, the most natural thing for Quintessence is to have a time varying equation of state, however one that satisfies 
the condition dlOt . 

The two models reviewed in this section, though very different, are inside the framework of General Relativity (GR). 
A further dramatical departure is to modify the laws of gravity and to consider DE models outside the framework of 
GR. We will consider in details the example of scalar-tensor DE models. 



SCALAR-TENSOR DARK ENERGY MODELS 



Full Reconstruction 

As we have seen, in quintessence models the equation of state parameter i 
inequality must be satisfied 
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must satisfy w^>—\. The following 
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Note that inequality dl 2l i applies only to spatially flat universes, a full analysis should relax the flatness prior 1 13']. It is 
not clear from the existing data whether fl!2i is satisfied and actually the analysis of the most recent SN data supports 
a varying equation of state which is of the Phantom type, i.e. with w < — 1, on very small redshifts < z < 0.5. 
If confirmed, a striking consequence is that Quintessence models are ruled out. It is therefore important to consider 
a more general class of models, like scalar-tensor (ST) models, where the inequality d!2i is no longer compulsory. 
Further ST theories are interesting to consider as they arise naturally from more fundamental theories like M-theory. 

We now review the reconstruction program to scalar-tensor DE models |9]. We consider the following Lagrangian 
density in the Jordan frame (JF) 



L=-[F(*)R-gV v d ll *d v Q)-U(4>)+L l , 



(13) 



where L m describes dustlike matter and F(<&) > 0. The Lagrangian as it is written in d 1 3I > can describe consistently 
models with a positive Brans-Dicke parameter C0bd — F / (dF / d<P) 2 > 0. If (Obd < 0, we must use the parametrization 
with Z = — 1 where Z is the function in front of the kinetic term, the Lagrangian dl 31 corresponds to the Z = 1 
parametrization. We will write all our equations using dl 3i . i.e. in the JF using the Z = 1 parametrization. We do 
not introduce any direct coupling between <t> and L m so that, in particular, fundamental constants do not change with 
time. We can proceed in a way analogous to the Quintessence case, the main difference being that we have now 
to reconstruct two unknown functions instead of one. Again, specializing to a flat FRW universe, the corresponding 
(modified) Friedmann equations read 
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As was the case in GR, the evolution equation for the scalar field <t> is contained in the two Friedmann equations above. 
These background equations d 1 4-1 1 5 i can be combined to give the following master equation for F(z): 
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where a prime denotes a derivative with respect to z. 

In this theory, the effective value of Newton's gravitational constant Gn is given by 
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It is natural to use its present value Gfjfi = 5^ in the definition of the critical density p cr jt- However Gn does not have 
the same physical meaning as in GR, it is no longer the coupling constant for the gravitational attraction between two 



point masses. For a massless dilaton, the effective gravitational constant between two test masses is given by 



(19) 



%kF \2F + 3{dF/d<P) 2 
In our case, the dilaton is massive but (II 81 will still hold for physical scales R such that 
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The most recent solar system measurements 1 14] imply very stringent constraints on the Brans-Dicke parameter 
today 

fd<t>\ 2 a 

(Obd,o = F Q i^—j > 4 x 10 4 . (20) 

As a consequence, G#o ar, d Geftfi coincide with better than 1.25 x 10~ 5 accuracy. On the other hand, the difference 
between Gn and G e ff could be larger at higher redshifts. Our interest for scalar-tensor theories of gravity is in the 
context of DE models. So our theory should satisfy the following requirements as any realistic DE model. First, the 
DE term should dominate today the energy density of the universe and satisfy 

£i D£ ,o ~ 0.7 ~ 2£2 m ,o . (21) 

If our model describes a universe whose expansion is presently accelerated, then it must satisfy 

U Q > (pm + 2* 2 + 3F + 3HF) (22) 

Finally it is important that DE remains essentially unclustered at scales up to R ~ 10/j 1 ( 1 + z) _1 Mpc, though as we 
will see some limited clustering will arise on these scales. To achieve this, it is sufficient to assume that the inequality 
dl9l is satisfied for all scales of interest. 

Like in GR, we start with the determination of H(z) from Di(z) using 
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However, we need to recover the two functions F(z) and U(z), so the substitution of H(z) in (I16> is no longer sufficient. 
For a complete reconstruction we must use a new equation based on independent observations. It is provided by 8,„(z) 
data, measurement of the linear dustlike matter density fluctuations. We can expect in the near future accurate Di(z) 
and 8 m (z) data. 

We consider the perturbation equations in the longitudinal gauge 

ds 2 = -(l+2((>)dt 2 + a 2 (l-2\i/)dx 2 , (24) 

(see for details). The idea is that, in the short wavelength limit, the leading terms are either those containing k 2 , or 
those with 8 m . Then the following equation is obtained 

. dF F dF/d<P 

Ws »- 2 ^» s -*fWw- (25) 

Hence, unlike in GR, in ST gravity the dilaton remains partly clustered even for arbitrarily small scales, however this 
clustering is small because <x>bd is large. In the same short wavelength limit, Poisson's equation has the same form as 
in GR, with the important difference that Newton's constant Gn is replaced by G e ff, defined in (II 81 above. Hence, the 
equation for the evolution of dustlike matter linear density perturbations finally leads to 

H 2 5: 8' m * hi +z)H 2 ^ £l m , 8 m . (26) 
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So, we see here a second difference in the physical meaning of G e ff and Gn: while it is Gn that appears in the equation 
for the evolution of the perturbations in GR, in ST models this role is played by G e g. 



Let us now sketch briefly the reconstruction itself. Extracting H(z) (through Dt,(z)) and 8 m (z) from observations 
with sufficient accuracy, we first reconstruct G e ff(z)/G^o analytically. Since, as follows from Eq.j20K the quantities 
Geff.o and Gn.q coincide with better than 0.00125% accuracy, Eq.(l26> taken at z = gives also the value of Q. m q with 
the same accuracy. Thus, in principle, no independent measurement of Q m o is required, it follows from ( I26> taken at 

z = 

We get an equation G e ff(z) = p(z), where p{z) is a given function that can be determined solely using observational 
data, which can be transformed into a nonlinear second order differential equation for F{z) using the background 
equation 
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Hence F(z) can be determined by solving the equation G e g- (z) = p{z) after we supply the initial conditions Fq = % k q nq 

and Fq. Actually F Q ' must be very close to zero due to the solar system constraint i20\ . Once F(z) is found, it can be 
substituted into equation ( 1161 to yield the potential U(z) in function of redshift. Then, using Eq. ( I27> . <t>(z) is found by 
simple integration which, after inverting this relation, gives us z — z( < t > — $q). Finally, both unknown functions F(<P) 
and U (O) are completely fixed as functions of <t> — <t>o- Of course this reconstruction can only be implemented in the 
range probed by the data corresponding to z < 2. 



Partial reconstruction 

As one does not expect to have in the very near future data referring to perturbations that are as accurate as the 
distance-luminosity data, it is interesting to try to extract as much information as possible using only £>l(z) data. If we 
make assumptions on either F(z) or U(z), or if we assume some functional relation between both, it is again possible 
to reconstruct the theory using only Dz(z) data. This is what we mean by reconstruction of constrained models, or 
partial reconstruction: we reconstruct models where by assumption there is effectively only one unknown function. 
Several cases have been considered in 1 10] and powerful constraints can be obtained. In particular, it is interesting that 
these constraints can go beyond solar-system constraints just because they go back in time and probe the cosmological 
evolution. To see how powerful constraints can be obtained, we consider the interesting question whether it is possible 
to have a vanishing potential U while the expansion of the universe still satisfies 

H 2 (z) =// 2 [0.3 (l+z) 3 +0.7] . (28) 

Of course in the framework of GR an expression like J2 81 lends itself to the straightforward interpretation of a flat A 
dominated universe with the cosmological parameters 

Ovo w 0.7 , a mfi « 0.3 . (29) 

The dynamical laws of GR are encoded in the definition of these cosmological parameters through the Friedmann 
equations. Now, we assume we have the same kinematics (expansion) but with modified dynamical laws as we deal 
with a modified gravity theory. When U = equation d!6t can be solved numerically and one can show that F(z) will 
vanish at most at Zmax ~ 0.66 if F(z) evolves according to dl6l . Therefore, the predicted vanishing of F(z) is well 
inside the range for which H(z) is determined with sufficient accuracy. Hence, for H(z) obeying d28t i up to z ~ 0.66, 
these models with vanishing potential U are excluded. 

If observations force us to consider models that can have DE of the phantom type then quintessence models, 
i;e. models inside General Relativity with a minimally coupled scalar (quintessence) field, are ruled out and is led 
to consider models outside GR.like scalar-tensor DE models that we have been discussing here. Many problems 
appearing in these models are representative of all DE models outside GR. Clearly discrimination between all the 
models will accurate data probing the background and the inhomogeneities: Supernovae data, Cosmic Microwave 
Background data, galaxy surveys, weak lensing data, etc. (see e.g. 11511 '). 

Recently, investigation of f(R) modified gravity theories led to the surprising result fl6ll that the cosmological 
history for a large class of such models is incompatible with observations because of the disappearance of the usual 

2 1 

matter-dominated stage with a^t^ which is replaced by a stage with a oc t 2 . So in this case, cosmological background 
constraints come from large z behaviour and not just from SNIa data on small redshifts z < 2! 
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